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CHAPTER  I 


INTRODUCTION 


This  report  presents  a  moment  method  (MM)  [1]  [2'  solution  to  two-dimen¬ 
sional  (2D)  transverse  magnetic  (TM)  scattering  by  a  variable  sheet  impedance  in 
a  multilayered  slab.  An  integral  equation  is  derived  for  the  sheet  impedance  in  the 
multilayered  slab.  This  integral  equation  is  solved  for  the  surface  currents  flowing 
on  the  sheet  impedance  by  employing  a  spectral  domain  MM/Green’s  function 
solution  [3].  This  solution  is  preferred  because  the  Green’s  function  acrounts  for  the 
presence  of  the  multilayered  slab.  Also,  the  fields  in  each  layer  of  the  multilayered 
slab  can  be  easily  expanded  as  plane  wave  spectra,  thus  the  MM/Green’s  function 
solution  is  carried  out  in  the  spectral  domain.  A  user  oriented  computer  code  was 
written  to  implement  this  solution  and  numerical  results  are  presented. 

A  sheet  impedance  is  a  model  for  an  electrically  thin  dielectric  layer.  The  sheet 
impedances  considered  here  can  be  variable,  but  they  must  be  non-constant  over 
a  single  finite  range.  The  multilayered  slab  is  planar  and  extends  infinitely  in  all 
directions.  The  slab  layers  are  of  dielectric/ferrite  materials  which  may  be  either 
lossy  or  lossless.  The  geometry  of  the  sheet  impedance  and  the  multilayered  slab 
is  shown  in  Figure  1  in  Chapter  II.  The  problem  is  to  compute  the  2D  scattering 
by  this  geometry  from  a  TM  incident  plane  wave. 

The  MM  solution  presented  in  this  report  is  useful  in  that  it  allows  for  the 
analysis  of  scattering  from  isolated  scattering  points  and  variable  sheet  imped- 


antes.  For  example,  scattering  from  a  single  sheet  impedance  discontinuity  can  lie 
computed,  as  well  as  scattering  from  a  tapered  sheet  impedance  where  the  lapei 
occurs  over  a  finite  width.  Also,  scattering  from  a  finite  width  sheet  impedance 
variation  can  he  computed.  Finally,  the  sheet  impedance  may  he  located  in  a 
mult  i  layered  si  ah,  on  a  half  plane,  or  ent  i  rely  in  free  space,  and  t  he  scat  t  eri  ng  from 
these  configurations  can  he  computed. 

The  MM  solution  is  begun  by  considering  the  sheet  impedance  and  multi 
layered  slab  to  be  illuminated  by  a  I'M  plane  wave  incident  from  the  angle  r/>o . 
The  variable  sheet  impedance  is  replaced  by  the  parallel  combination  of  two  sheet 
impedances.  One  sheet  impedance  will  be  constant  while  the  other  sheet  imped¬ 
ance  will  be  variable  and/or  constant  and  will  produce  the  scattered  field.  Next, 
the  equivalence  theorem  is  used  to  replace  the  “scattering"  sheet  impedance  by 
an  unknown  surface  current  J,.  The  condition  that  the  total  field  is  the  sum  of 
the  incident  field  and  the  scattered  field  yields  the  integral  equation  for  Js.  The 
unknown  surface  current  Js  is  expanded  in  terms  of  N  known  basis  functions  of 
unknown  strengths.  The  /V  coefficients  in  this  expansion  are  then  determined  by 
the  moment  method.  The  scattered  field  is  the  sum  of  the  individual  fields  of  the 
N  basis  functions  multiplied  by  their  respective  strengths. 

Some  previous  work  has  been  published  on  the  impedance  sheet  approxima¬ 
tion  and  resistive  and  impedance  sheet  scattering.  Harrington  and  Mautz  |l|  pre 
seated  the  impedance  sheet  approximation  and  applied  it  in  a  MM  solution  for  TM 
scattering  by  a  thin  dielectric  shell.  Senior  [fij  computed  edge-on  I'M  backscatter 
ing  from  a  uniform  resistive  sheet  half  plane  in  free  space.  Senior  jfil  also  analyzed 
back  scattering  from  finite  widt  h  resistive  st  rips  with  attention  focused  on  I'M  scat 
tering  near  grazing  incidence.  Again,  his  analysis  concerned  strips  with  uniform 
resistance.  Richmond  |7|  presented  a  MM  solution  for  I'M  scattering  by  a  finite 


width  electrically  thin  dielectric  strip.  II is  solution  used  basis  functions  incorpo¬ 
rating  known  physical  properties  of  the  problem  geometry  into  the  MM  solution. 
The  MM  solution  presented  in  this  report  makes  use  of  similar  basis  functions. 
Newman  [8]  [9]  used  a  MM /Green's  function  solution  to  solve  for  scattering  bv  a 
dielectric/ferrite  cylinder  in  the  presence  of  a  perfectly  conducting  half  plane.  This 
solution  can  be  used  to  model  a  l’ EC  half  plane  with  a  variable  impedance  taper 
region  at  its  edge.  Newman  and  Blanchard  [10]  solved  for  I'M  scattering  by  an 
impedance  sheet  extension  of  a  parabolic  cylinder  using  a  MM/Greens  function 
solution. 

Asymptotic  work  employing  the  uniform  theory  of  diffraction  (lTTD)  (I  l]  has 
also  been  done  in  the  analysis  of  thin  dielectric/ ferrite  slabs.  Rojas  and  Pathak 
[12]  [13]  analyzed  diffraction  by  dielectric/ferrite  half  planes  and  strips.  Rojas  [1  1] 
also  solved  for  scattering  by  an  impedance  discontinuity  in  a  planar  surface.  Ly 
(15f  presented  a  UTD  solution  for  diffraction  by  junction  edges  formed  between 
different  electrically  thin  material  slabs.  For  simple  cases,  the  MM  solution  pre¬ 
sented  here  was  compared  against  the  solution  given  in  [15].  The  agreement  was 
always  excellent. 

The  format  of  this  report  is  as  follows.  The  basic  geometry  of  the  problem 
is  given  in  Chapter  II.  Chapter  II  also  derives  the  integral  equation  and  outlines 
the  MM  solution  for  the  sheet,  impedance  in  the  multilayered  slab.  Chapter  III 
presents  results  of  TM  plane  wave  scattering  by  several  different  sheet  impedances 
and  multilayered  slab  geometries.  When  possible,  these  results  were  compared  with 
measured  or  previously  calculated  results.  Chapter  IV  briefly  describes  a  computer 
code  written  to  implement  this  MM  solution  and  is  intended  as  a  user's  manual 
for  the  code.  Finally,  Chapter  V  gives  a  summary  of  this  report,  and  presents  ideas 
for  further  study. 
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CHAPTER  II 


THEORY 

2.1  Introduction 

I  his  chapter  develops  the  integral  equation  and  moment  method  solution  in 
the  problem  of  2 It  transverse  magnetic  (I'M)  plane  wave  scattering  by  a  variable 
sheet  impedance  in  a  planar  multilayered  slab.  I'he  sheet  impedance  must  be  non- 
constant  over  a  finite  range  for  the  MM  sob-lion  presented  here.  I'he  slab  geometry 
is  illustrated  in  figure  I  and  a  sheet,  impedance  model  is  shown  in  Figure  2. 

The  solution  is  obtained  by  first  expressing  the  variable  sheet  impedance  as 
the  parallel  combination  of  two  sheet  impedances.  One  sheet,  impedance  will  be 
constant  while  the  other  sheet  impedance  will  be  infinite,  constant  or  variable 
over  different  ranges  and  produces  the  scattered  field.  The  equivalence  theorem  is 
used  to  replace  the  scattering  sheet  impedance  by  an  unknown  surface  current  ,/s. 
Next,  the  integral  equation  for  ./y  is  obtained  by  enforcing  the  surface  equivalence 
theorem  on  the  sheet  impedance  surface.  I  bis  integral  equation  is  solved  using  a 
spectral  domain  moment  method  ( M  M  )/Creen?s  function  solution  which  solves  for 
the  surface  current  ,/.s.  Once  is  known  then  the  scattered  fields  are  computed. 

2.2  Problem  Geometry  ami  the  Integral  Equation 

I  his  sect  ion  present  s  the  development  of  an  in  I  egral  equa  t  ion  for  the  2 1 )  t  ran< 
verse  magnetic  (  IM)  scattering  by  a  variable  sheet  impedance,  denoted  /s(.r). 


I 


in  a  plane  multilayered  slab.  As  illustrated  in  Figure  1,  the  slab  lias  h  layers 
with  permeability  and  permittivity  (//£,<&)  and  thickness  7).;  k  1.2,....  /v  .  l'he 
material  parameters  of  the  free  space  region  above  the  multilayered  slab  are  de¬ 
noted  (p Q .  r o )  and  the  homog^^ous  region  below  the  multilayered  slab  are  denoted 
(/*A'4  l>cA't  1  )•  AM  fields  and  currents  are  two  dimensional  (they  are  independent 
of  y )  and  are  time  harmonic  with  the  cJU><  time  variation  suppressed.  In  addition, 
all  electric  fields  and  surface  currents  contain  only  a  y  -  directed  component  so 
vector  notation  shall  not  be  used  and  Ibis  polarization  is  implicit. 

A  sheet  impedance  is  a  model  for  an  electrically  thin  dielectric  slab.  For 
example,  Figure  2(a)  shows  a  dielectric  slab  of  thickness  T  and  with  material 
parameters  (/q),f).  The  wavenumber  in  the  slab  is  k  =  If  \k\T  1, 

then  the  slab  is  sufficiently  thin  that,  the  electric  field  is  essentially  constant  with 
respect  to  s.  In  this  case,  as  illustrated  in  Figure  ‘2(b),  the  thin  dielectric  slab  can 
be  replaced  by  the  zero  thickness  sheet  impedance  (I0j 


ZS 


7^0 


<o)T 


(2.1) 


Although  Figure  2(a)  shows  a  homogeneous  dielectric  slab  of  uniform  thickness, 
Equation  (2.1)  also  applies  to  a  slab  in  which  T  and/or  e  are  functions  of  x. 

The  sheet  impedances  considered  lure  have  the  property  that  for  some  ranges 
of  x ,  Z$(x)  ~  Z$ o  is  constant,  while  for  other  ranges  of  x,  Zc;(x )  is  an  arbitrary 
function  of  x.  Furthermore,  the  range  over  which  Z$(x)  is  non-constant  (Z50  or 
otherwise)  must  be  finite  t  .>  allow  for  a  feasible  sub-domain  basis  function  expan¬ 
sion  in  this  region.  Some  typical  Zy(,r)  are  shown  in  Figure  3.  Figure  3(a)  shows 
the  most  general  sheet  impedance  with  an  arbitrary  variation  for  0  <  x  <  L  and 


layers  (k+3)  -  (K) 

- j 

homogeneous  (m 

K+l ^  6  K+l) 


Figure  I:  Mull  Havered  sla.l>  and  sheet  impedance  geomet  ry. 


Etan  =  total  tangential  electric  field 

(b) 


Figure  2:  Sheet  impedance  approximation  through  the  equivalence  theorem. 


Z.s(.r) 


which  is  constant  at  Z$ q  or  Z  j  outside  this  range.  For  example,  in  Figure  .'{(h) 

' 

Zsu  if  x  <  0 

Zs(*)  {  Z,  iff)<;r  <  /,  (2.1’) 

Zs„  if  x  >  L. 

The  Z^(x)  in  Figures  .'{(c)  and  3(d)  differ  from  Zy( )  over  a  semi  infinite  range. 
However,  Z$ ( x )  is  11011-constanf  over  a  finite  range  of  ,r,  as  required.  In  particular, 
for  Figure  3(c) 

( 

Z s"o  if  x  <  0 

(2.3) 

Z,\  if  x  >  0, 

\ 

and  for  Figure  3(d) 

Z.vo  if  x  <  0 

^S'(J-)  Z i\x)  if  0  <  r  <  L  (2,1) 

Z  .|  if  a-  >  L 

where  Zj'(x)  is  an  impedance  taper  function  providing  a  smooth  transition  from 
ZSl]  to  Z  4  over  the  range  0  <  x  L. 

A  simplified  sketch  of  the  multilayered  slab  containing  Zy(a’)  is  shown  in 
Figure  -1(a).  Newman  showed  that  the  equivalent  sheet  impedance  of  a  thin  multi¬ 
layered  dielectric  slab  is  given  by  the  parallel  combination  of  the  sheet  impedances 
of  the  individual  layers  (10).  Thus,  the  sheet  impedance  Zy(  x )  can  be  represented 
as  the  constant  sheet  impedance  Z<,u  in  parallel  with  another  sheet  impedance,  de 
noted  Z$p(x).  This  equivalent  geometry  is  shown  in  Figure  -1(b).  Z^(ai)  is  given 
by 

^.sd%/’(;r) 


^S{X)  y  |  7  /  \ 

Z‘S\  I  1  /'Sr{x) 

from  which  it  is  found  that  Zy />(./■ )  is  given  by 


Zsp(x) 


ll  ■  7s(x) 


(2.r>) 


(2.h) 


8 


Figure  I :  hqui valent  geometries  used  iu  deriving  the  integral  equation. 


11) 


The  scattered  field  is  produced  by  the  parailel  sheet  impedance  Zgp(x),  which 
is  finite  only  where  Zg(x)  ^  ZgQ.  As  illustrated  in  Figure  4(c)  the  volume  equiv¬ 
alence  theorem  can  be  used  to  replace  ZgP(x)  by  the  surface  current  [10] 


Js(*) 


El 

ZSp(x)' 


(2.7) 


where  Et  is  the  total  electric  field  on  the  surface  of  Zgp(x).  Note  that  J s{x )  is 
non-zero  only  where  Zg( x)  /  Zg q,  and  that  Jg{x)  radiates  the  scattered  field  Es 
in  the  presence  of  the  multilayered  slab  of  Figure  4(c).  The  term  “multilayered 
slab”  refers  to  the  K  material  layers  and  the  sheet  impedance  ZgQ,  but  not  Zgp(x). 
The  total  electric  field  is  the  sum  of  the  incident  plus  the  scattered  field,  i.e., 


El  =  EJ  +  Es 

where  El  is  the  field  of  the  incident  plane  wave  in  the  presence  of  the  multilayered 
slab  of  Figure  4(c).  Combining  Equations  (2.7)  and  (2.8)  yields 


(2.8) 


-  Es  +  ZSP(x)Jg(x)  =  E\ 


(2.9) 


This  is  the  basic  integral  equation  for  J g(x)  because  the  electric  field  Es  can  be 
expressed  as 

Es  =  jCsy(x'\x,z)Jg(x')dx'  (2.10) 

w'here  Gy(x  ]x,z)  is  the  y  component  of  the  dyadic  Green’s  function  for  the  mul¬ 
tilayered  slab  containing  ZgQ.  The  y  component  of  the  dyadic  Green’s  function  is 
chosen  because  the  surface  current  is  y  -  directed.  However,  the  dyadic  Green’s 
function  will  not  be  used  in  this  analysis  but  is  used  only  to  show  the  form  of  the 
integral  equation.  Equation  (2.9)  will  be  solved  for  Jg(x)  in  the  next  section  using 
the  MM. 
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2.3  Moment  Method  Solution 


This  section  presents  the  development  of  the  moment  method  solution  to  the 
integral  equation  developed  in  the  previous  section.  The  unknown  surface  current 
is  approximated  by  a  sum  of  /V  known  basis  functions  of  unknown  strengths.  1'he 
choice  of  the  basis  functions  is  discussed  in  Section  2.8.  By  taking  an  inner  product 
of  Equation  (2.9)  with  a  set  of  N  weighting  functions,  a  matrix  equation  is  formed. 
The  weighting  functions  will  be  chosen  identical  to  the  basis  functions  making  this 
a  Galerkin  solution  [ J ] .  Finally,  general  expressions  for  the  matrix  elements  are 
given.  Note  that  all  electric  fields  and  surface  currents  are  y  -  directed  so  vector 
notation  is  not  used. 

The  moment  method  solution  is  begun  by  approximating  the  unknown  surface 
current  Jsix)  as  a  sum  of  N  known  basis  functions  as  follows: 

N 

JS(x)  Y  ln  Jn  (2.11) 

n  ~ .  I 

where  the  In  are  the  N  unknown  strengths  of  each  basis  function  and  the  Jn  are  the 
N  known  basis  function  expansion  modes.  Substituting  the  approximate  surface 
current  Js{x)  Equation  (2.1  I)  into  Equation  (2.9)  it  is  obtained  that 

N  N 

Y  In  En  t  V  2.S'M*)  &  (2.12) 

71  I  11  I 


where  En  is  the  electric  field  of  the  n1^  basis  function  ,/„  radiating  in  the  presence 
of  the  multilayered  slab  (containing  /f.so).  Bv  taking  the  product  of  Equation 
(2.12)  with  each  weighting  function  for  m  1,2,...,  /V,  and  integrating  then 
over  the  region  of  the  weighting  function,  it  is  obtained  that 


yv  n 

Y,  In  /  EnJmdr  +  Y  1»  /  • ln-l,„Zsy(.r)dx  /  K1  Jtndr.  (2.12) 

.  Jm  i  >hn  Jm 
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The  limit  m  on  the  integrals  in  the  above  equation  denotes  that  the  range  of 
integration  is  over  the  region  where  thJm  is  non-zero.  Equation  (2.13)  represents 
a  set  of  N  equations  with  N  unknowns.  This  set  of  equations  can  be  written 
compactly  in  matrix  form  as 

|2  +  AZ||f|  =  [l’|  (2.1.1) 

where  [Z  +  AZj  is  the  N  x  N  impedance  matrix,  [/]  is  the  length  N  vector  of  un¬ 
known  strengths,  and  [H]  is  the  length  N  voltage  excitation  vector.  The  following 
equations  define  the  elements  of  the  impedance  matrix  and  the  voltage  vector: 

Zmn  ~  f  EnJrndx  (2.15) 

Jm 

—  [  Jn  J m  Z gp(x)  dx  (2.1G) 

Jm 

I'm  —  f  El  Jmdx.  (2.1  7) 

Jm 

Reciprocity  can  be  applied  to  Equation  (2.17)  to  obtain  an  expression  for  the 
voltage  vector  elements  in  terms  of  Em,  the  electric  field  of  the  mth  weighting 
function  Jm.  The  resulting  reciprocal  expression  is 

Vm  —  J  Emji  dv  (2.J  8) 

where  Jl  is  the  impressed  current  that  radiates  the  incident  electric  field  E\  The 
integral  is  over  the  volume  of  the  impressed  current. 

The  Zmn  impedance  matrix  contributions  are  distinctly  different  from  the 
A Zmn  contributions.  Equation  (2.15)  contains  the  electric  field  of  the  basis  func¬ 
tions  radiating  in  the  presence  of  the  multilayered  slab.  However,  Equation  (2.16) 
does  not  contain  any  electric  fields.  The  Zmn  terms  result  from  the  field  of  the 
basis  functions  in  Equation  (2.9),  whereas  the  A Zmn  terms  result  directly  from 
the  basis  functions.  As  a  result,  the  Zmn  terms  require  further  evaluation,  but  the 


AZjjui  terms  can  be  evaluated  directly  from  Equation  (2.lfi)  for  known  basis  and 
weighting  functions  and  surface  impedance. 

Note  that  both  Equations  (2.15)  and  (2.18)  contain  the  electric  fields  of  I  la- 
basis  functions  and  weighting  functions  radiating  in  the  presence  of  the  multilay¬ 
ered  slab.  The  next  section  develops  expressions  for  these  electric  fields  and  uses 
them  to  evaluate  the  impedance  matrix  elements. 

2.4  Evaluation  of  the  Impedance  Matrix  Elements 

This  section  obtains  an  exact  integral  expression  tor  the  impedance  matrix 
elements  from  the  general  expression  given  by  Equation  (2.15)  in  the  proceeding 
section.  An  expression  for  the  total  electric  field  of  the  basis  function  radiating 
in  the  presence  of  the  multilayered  slab  will  be  developed  using  the  plane  wave 
expansion  method  |l(i|.  This  expression  will  be  used  in  the  evaluation  of  tin- 
matrix  elements  of  Equations  (2.15)  and  (2.18).  This  section  evaluates  Equal  ion 
(2.15)  in  the  spectral  domain. 

A  separate  plane  wave  expansion  for  the  fields  is  used  in  each  layer  of  the 
slab,  in  the  free  space  region  above  the  slab,  and  in  the  homogeneous  region  below 
the  slab.  Thus,  for  a  A  layered  slab,  there  will  exist  A  (  2  separate  regions,  each 
with  its  own  plane  wave  expansion.  In  the  following  analysis,  the  notation  for  the 
subscript  and  superscript  k  is: 

•  k  0  implies  the  free  space  region  above  the  multilayered  slab, 

•  k  1,2,...,  A  implies  t  he  k^'  layer  of  the  multilayered  slab,  and 

•  k  All  implies  the  homogeneous  region  below  the  multilayered  slab. 


II 


The  basis  functions  Jn  are  assumed  to  be  Fourier  transformable  with  the 


Fourier  transform  pair  defined  as 


J,j  -  [  .7,i  e±Jf3x  dx 

J  -OO 


+  °°  j 7  jd* 

-  oo 


'  n 


dfl. 


(2.19) 

(2.20) 


are  referred  to  as  the  “plus”  and  “minus”  Fourier  transforms  of  Jn. 


The  surface  current  lies  on  the  interface  between  two  regions  so  each  region 


can  be  considered  source  free.  Thus,  the  electric  and  magnetic  fields  in  each  region 
must  satisfy  the  source  free  vector  wave  equation 


(v2  f  4) 


(2.21) 


where  k j.  =  's  ^,e  wave  number  of  the  kth  region.  The  total  y  -  po¬ 

larized  electric  field  in  region  k,  produced  by  Jn  radiating  in  the  presence  of  the 
multilayered  slab,  is  expressed  as  a  continuous  spectrum  of  plane  waves,  i.e., 


E*{X,  z)  -  r  °°  [  A* (/?)  +  B*(/9)  e^z  1  e^x  d/3  (2.22) 

where  A^(0)  and  B^(0)  are  spectral  weighting  functions  to  be  determined,  and 
is  a  propagation  constant  for  waves  travelling  in  the  z  direction.  Applying 
Equation  (2.21)  to  Equation  (2.22)  yields  the  separation  equation 


7 k  =  \/^2  -  kl  ;  Re(7 jfc)  >  0,Ini(7jfc)  >  0.  (2.23) 


This  relates  the  propagation  constant  7*.  to  the  plane  wave  spectrum  parameter 
0.  The  conditions  on  7^  insure  that  the  wave  decays  in  the  direction  it  travels. 
In  this  manner,  the  term  of  Equation  (2.22)  associated  with  AXI(0)  is  an  upward 
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travelling  wave  and  the  term  associated  w  ith  #*(/?)  is  a  downward  travelling  wave. 
Expressions  for  the  magnetic  field  are  obtained  from  Maxwell’s  equation 

V  x  E  (2.21) 

Applying  Equation  (2.21)  to  Expiation  (2.22)  the  following  magnetic  field  is  ob- 
I  ai  ned : 


|  _ 

??  i 

pt- 

l-lfilXIr  ’*  = 

nt(ri)<lkT  |  >J-h'  ,/n 

( 2.2b ) 

—  r~,i  i 

“7'  it-  '  ™ 

1  1 

(  2 . 2  <  i ) 

Next,  tlic  functions  .l^(.rf)  and  /t^’(if)  must  he  determined.  These  functions 
are  determined  by  enforcing  tin'  following  conditions  on  the  (ields: 

1.  the  radiation  condition  as  r  •  x  . 

2.  continuity  of  the  tangential  electric  held  across  each  interface. 

!{.  continuity  of  the  tangential  magnetic  field  across  each  interface  except  the 
:  0  interface,  and 

I.  discontinuity  of  the  tangential  magnetic  field,  by  the  total  surface  current 
produced  by  across  the  r  U  interface,  i.e., 

|nf,  H* 1  1  j  j'„  (2.27) 

where  b  and  b  I  I  represent  the  regions  immediately  above  and  below  the 
z  (1  interface  respectively,  and  .l(,  is  the  total  surface  current  produced  by 


Ifi 


jJ  contains  only  a  y  -  directed  component  which  at  the  2-0  interface  can  be 


written  as 


%(2') 


Jn 


En{x,0) 

^50 


(2.28) 


In  Equation  (2.28)  E% +  1  could  have  been  used  instead  of  E *  since  the  tangential 
electric  field  is  continuous  across  the  z  -  0  interface. 

Enforcing  the  above  conditions  on  the  fields  of  Equations  (2.22)  and  (2.25), 
and  substituting  Equations  (2.20)  and  (2.22)  into  Equation  (2.28),  the  electric  field 
E k  is  obtained  in  terms  of  In  Appendix  A  it  is  shown  that  the  electric  field  at 
the  2  =  0  interface  is  given  by 


*£(*,«)  -  ~  F<'a>  ■'» e>01  &  (2-2!)) 

where  F(/3)  is  a  spectral  function  specific  to  the  multilayered  slab  geometry.  E(/3) 
can  be  viewed  as  a  spectral  domain  Green’s  function.  See  Appendix  A  for  the 
determination  of  F{(3)  for  several  simple  geometries,  such  as  the  sheet  impedance 
located  in  a  K  =2  layered  slab,  on  a  material  half  space,  and  entirely  in  free  space. 

The  impedance  matrix  elements  are  obtained  by  substituting  Equation  (2.29) 
into  Equation  (2.15)  to  obtain 

Zmn  I  F(/3)  j-  erf*  Jrn  d(3  dx.  (2.30) 

2  7T  Jm  J  —  oo 

Reversing  the  order  of  integrations,  and  employing  the  Fourier  transform  relation¬ 
ship  of  Equation  (2.19),  it  is  finally  obtained  that 

Zmn  ^  ^  F(0)  j+  J-  d0.  (2.31  ) 

This  is  the  so-called  spectral  integral  formula  for  the  impedance  matrix  elements. 
This  integral  expression  is  preferable  to  using  Equation  (2.10)  in  Equation  (2.15) 
because  there  is  only  one  integration.  Furthermore,  the  dyadic  Green’s  function 
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for  the  mull  ilnvered  slab  need  no!  In*  del  ermined.  However,  1  lie  fiiuetioii  /'(d) 
mest  he  (letennined.  Note  dial  determining  (lie  impedance  matrix  elements  for 
different  multilayered  slab  geometries  is  trivial.  All  that  needs  to  he  done  is  to 
insert  the  /id)  corresponding  to  the  f'o'iiirlrv  of  interest  into  Equation  (2.21). 

I  Inis,  different  spectral  functions  /(d)  can  be  used  for  flilh’i'ent  geometries. 

2.5  Numerical  Evaluation  of  the  Impedance  Matrix 

The  impedance  matrix  elements  given  by  Equation  (2.21)  will  he  evaluated 
numerically.  There  are  two  problems  encountered  in  the  numerical  evaluation 
of  the  impedance  matrix.  The  first  problem  is  the  CIM1  time  required  for  the 
compulation  of  the  impedance  matrix.  The  second  problem  is  that  surface  wave 
poles  produce  singularities  in  the  integrand  of  the  spectral  integral  formula.  I  his 
section  discusses  these  problems. 

2.5.1  CPU  Time  Reduction 

The  spectral  integral  formula  for  the  impedance  matrix  elements  is  given  by 
Equation  (2.21).  Noting  the  form  of  the  equation,  the  re-evaluation  of  the  same 
quantities  can  be  eliminated  by  computing  all  the  elements  in  a  parallel  fashion. 
Thus,  the  entire  impedance  matrix  should  be  computed  in  just  one  numerical 
integration  with  a  matrix  multiplication  performed  inside  the  integral,  i.e.. 

Z  1  I  #12  •  •  ■  #  I  ,V 

Zu  |  /■>■>  •  •  ■  /  >  y 

Z-Y  |  Z  \  L>  •  ■  Z  y  v 

(2.22) 

Note  that  in  this  evaluation  of  I  lie  impedance  matrix,  the  quantities  /’(d). 
and  ./„  are  evaluated  only  once  at  each  value  of  d  employed  in  the  numerical 

IN 


integration.  This  results  in  an  immense  savings  in  CPU  time  required  for  the 
impedance  matrix  computation. 

2.5.2  Surface  Wave  Poles 


The  spectral  function  F(f3)  in  the  expression  for  Zrnn  of  Equation  (2.31)  may 
contain  poles  which  make  the  integrand  singular.  Physically,  these  poles  corre¬ 
spond  to  surface  waves  that  exist  on  the  multilayered  slab  containing  the  constant 
sheet  impedance  Z$q.  The  number  of  surface  wave  poles  can  increase  with  in¬ 
creasing  frequency.  Typically,  for  lossless  layers,  these  poles  lie  on  the  Re(/3)  axis 

between  kn  and  kn  -  the  maximum  wavenumber  of  a  layer  medium.  For  lossy 
u  Umax 

layers,  these  poles  will  move  into  the  second  and  fourth  quadrants,  but  for  only 
slightly  lossy  layers  they  will  be  close  to  the  Re(/?)  axis.  If  the  poles  are  close  to  or 
on  the  Re(/?)  axis,  they  will  distort  the  numerical  (3  integration  for  Zmn  unless  the 
integration  path  is  staggered  around  them  as  shown  in  Figure  5.  This  staggered 
contour  will  produce  an  accurate  result  since  the  staggered  contour  does  not  cross 
any  poles,  nor  does  it  pass  too  closely  to  their  singularities.  This  staggered  contour 
method  has  worked  well  for  A  given  by 


A- 


0.2 


if  IT  <  5A  tj 

~ymin 


(2.33) 


An  /IF  otherwise 

^  1Jmm ' 

where  II  is  the  width  of  the  pulse  basis  function  region  and  A £>  -  's  the  wavelength 
in  medium  ^£>mar-  The  contour  then  returns  to  the  Re(/3)  ax's  around  1.2 kj)maT. 


2.6  Evaluation  of  the  Voltage  Vector  Elements 


This  section  obtains  an  expression  for  the  voltage  vector  elements  from  the 
expression  given  in  Equation  (2.17)  for  the  case  of  a  uniform  incident  plane  wave. 
The  geometry  of  the  incident  wave  is  shown  in  Figure  4(a).  An  impressed  current 


19 


p  -  plane 


Im  (p) 


k-n 


m  ai 


t 


/ 


Staggered  path  for 
numerical  integration 


t  Re 


(P) 


A- 


o  i  i  k 

%  J1W 1 

\ 

\ 

poles 


figure  5:  Integration  contour  for  Z,,tn  in  jd  plane. 

source  t  hat  radiates  a  plane  wave  in  the  far  zone  in  free  space  will  be  defined  so  i  hat 
the  reciprocal  expression  of  Imitation  (2.18)  can  be  used.  Finally,  an  expression  for 
the  far  zone  field  of  a  basis  function  radiating  in  the  presence  of  the  multilayered 
slab  will  be  obtained.  This  far  zone  field  expression  will  be  used  in  the  evaluation 
of  F.quation  (2.18)  and  also  in  the  computation  of  the  far  zone  scattered  field.  I  he 
formulation  outlined  here  assumes  that  the  wave  is  incident  from  the  upper  half 
space  c  >  (J.  A  similar  approach  can  be  used  if  the  wave  is  incident  from  the  lower 
half  space. 

The  incident  plane  wave  illuminating  the  multilayered  slab  is 


yt  7  *"(|t  4  sin<*<|) 


(2.21) 


is  field  will  be  produced  by  the  impressed  current  line  source 

XttI'i 


“'l>  a 


:2.2a) 


located  at  the  point  (p,</>y)  In  the  limit  as  p  —*  oo.  It  should  be  noted  that  the 
strength  of  Jl  depends  upon  the  distance  p,  so  that  Er  has  unit  amplitude  near 
the  origin. 

Substituting  Equation  (2.35)  into  Equation  (2.18)  it  is  obtained  that 

I'm  =  --  Em  y/peik°p  (2. 30) 

W/X() 

where  Em  is  the  field  of  Jm  evaluated  at  (p,<^>(j),  the  location  of  the  impressed 
current  source  Jl .  The  method  of  stationary  phase  is  used  to  asymptotically 
evaluate  Em  (17)  [18]  in  Appendix  B.  Em  is  found  to  be  of  the  form 

P 

Em  —  P m(</>l))  p 

VP 

where  is  a  function  of  </>o  dependant  upon  the  multilayered  slab  geometry. 

See  Appendix  B  for  the  evaluation  of  Em  for  several  simple  geometries.  Substitut¬ 
ing  Equation  (2.37)  into  Equation  (2.36)  the  voltage  vector  elements  are  found  to 
be 

Vm-  -Fruit o)^5e"jf.  (2.38) 

This  completes  the  discussion  of  the  calculation  of  the  voltage  vector  elements. 
2.7  Computation  of  the  Scattered  Field 


Once  the  impedance  matrix  and  voltage  vector  are  known,  the  current  vector 
can  be  obtained  using  standard  matrix  algebra.  Then,  once  the  current  vector  is 
known,  the  far  zone  scattered  field  can  be  obtained  as  (he  sum  of  the  far  zone  field 
of  each  basis  function  Jn  multiplied  bv  its  respective  strength  /„.  Thus,  the  far 
zone  scattered  field  is  be  given  by 


IV 


N 

v  T  fV 

/  ,  lnLjn 
n—  1 


(2.39) 
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where  /„  is  I  lie  strength  of  ./„  and  JJ  is  the  far  zone  fit-id  of  ./„  radiat  ing  in  1  he 
presence  of  tlu-  mult, ilayered  slab.  Tin-  far  zone  fit-id  <>f  a  basis  fnnrtion  is  deter 
mint’d  asymptotically  in  Appendix  If  and  is  given  by  Equation  (lf..ri).  Substituting 
Equation  (If. 5)  into  Equation  (2.311),  the  scattered  field  is  found  to  be 

N  r  )hu> 

K*  Y.  Inl'nit)— T~ 

n  —  1  y/P 

where  P, ,((/>)  is  dependant  on  the  multilayered  slab  geomet  ry  and  is  also  determined 
in  Appendix  If  and  is  given  by  Equation  (lf.fi). 

2.8  Basis  and  Weighting  Functions 

This  section  discusses  the  basis  and  weighting  functions  used  in  the  MM  so¬ 
lution.  The  weighting  functions  are  chosen  identical  to  the  basis  functions  making 
this  a  (Jalerkin  solution.  It  will  be  explained  how  the  choice  of  basis  functions 
is  determined  by  the  nature  of  the  sheet  impedance  Zy(.r).  All  basis  functions 
are  Fourier  transformable,  as  required  by  the  spectral  integral  formula  for  the  im¬ 
pedance  matrix  elements.  The  basis  functions  and  their  Fourier  transforms  are 
presented  below. 

Recall  that  the  surface  current  J s-(a-)  is  non-zero  in  the  range  of  X  where 
Zg(x)  ^  Zyq.  Therefore,  Jg(x)  must  be  expanded  in  terms  of  basis  functions  in 
this  range.  Furthermore,  the  range  where  Z$(t)  /  Zyq  can  be  either  finite  or 
semi-infinite.  However,  the  range  where  Zy(r)  is  non  constant  must  be  finite  to 
allow  for  a  feasible  subsecl  ional  basis  function  expansion  in  this  range. 

f  igure  :f(b)  shows  a  case  where  Zs(.r)  /  Zy n  over  a  finite  range  of  .r.  In  this 
case,  J s{x)  >s  expanded  in  subsectional  basis  functions  over  this  range.  Pulse  func¬ 
tions  are  chosen  as  the  subsectional  basis  functions  in  this  MM  solution.  Figures 
3(c)  and  3(d)  show  cases  where  Z$(r)  /  Zcq  over  a  semi-infinite  range  of  x.  In 
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these  cases  it  would  he  impossible  to  express  Jsix)  as  an  expansion  of  subsectional 
basis  functions  over  the  entire  semi-infinite  range  because  N  oo.  lo  circumvent 
this  problem,  Js(x)  >s  expanded  in  terms  of  pulse  basis  functions  only  in  and  near 
the  range  where  Z$(x )  is  non-constant  or  discontinuous,  fn  addition,  Js(r)  is 
expanded  in  terms  of  physical  semi  infinite  domain  basis  functions  outside  of  this 
range.  A  physical  basis  function  is  a  basis  function  with  a  variation  dependant  on 
a  physical  characteristic  of  the  problem  geometry.  These  basis  functions  incorpo¬ 
rate  known  variations  into  the  MM  solution.  Thus,  some  insight  into  the  problem 
is  required  to  choose  this  type  of  basis  function.  For  example,  a  physical  optics 
(PO)  basis  function  is  chosen  to  account  for  the  current  variation  induced  by  the 
incident  plane  wave.  Also,  surface  wave  basis  functions  are  chosen  when  surface 
waves  can  be  supported  by  the  sheet  impedance.  These  are  the  two  types  of  phys¬ 
ical  basis  functions  used  in  this  MM  solution.  They  are  semi-infinite  domain  basis 
functions  because  they  extend  to  infinity  in  the  -  direction.  It  is  assumed  that 
the  current  far  enough  away  from  any  non-constant  variations  or  discontinuities 
of  Zg(x)  can  be  approximated  accurately  by  only  the  PO  and  surface  wave  basis 
functions  because  Z$(x)  is  constant  in  this  range. 

Figure  6  shows  the  choice  and  placement  of  basis  functions  needed  for  a  proper 
and  complete  expansion  of  J g(x)  for  three  typical  variations  of  Zg(x).  It  should  be 
noted  that  in  Figures  6(b)  arid  6(c)  surface  wave  basis  functions  should  be  included 
only  if  the  multilayered  slab  geometry  can  support  surface  waves.  However,  surface 
wave  basis  functions  need  not  be  included  (even  if  the  geometry  can  support  them) 
if  the  surface  wave  decays  to  insignificant  strength  in  the  pulse  basis  function 
region.  If  the  slab  layer  is  lossy  enough,  or  the  surface  wave  propagation  constant 
provides  rapid  enough  attenuation,  then  this  will  be  the  case.  In  these  cases,  the 
current  variation  near  the  non-constant  or  discontinuous  sheet  impedance,  caused 
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by  t  lie  attenuated  surface  wave,  will  lie  approximated  well  enough  by  I  lie  pulse 
basis  functions  alone. 


If  1  lie  range  where  Zy(,r)  /  /su  it*  finite,  then  pulse  basis  functions  are  us<u| 
over  the  entire  range.  II  /<,-( j)  /  Zyu  over  a  semi  infinite  range,  then  pulse 
basis  functions  are  used  in  and  mar  the  range  tv  here  Zsf.;)  is  inai  constant  <<r 
discontinuous,  i.e..  they  are  used  over  a  range  starting  where  Zy(.r)  first  deviates 
from  Zy(),  and  ending  far  enough  beyond  where  has  become  constant  once 

again  that  Js(x)  tan  ^H'  approximated  by  only  the  PO  and  surface'  wave  basis 
functions.  In  either  case,  the  range  where  ■/<,•(  a1)  is  expanded  in  terms  of  pulse 
basis  functions  is  finite.  This  finite  range  is  divided  into  N  segments  and  a  pulse 
basis  function  is  placed  on  each  segment.  I'lie  pulse  basis  functions  are  of  the  form 


loss  vanishes,  then  Jn  will  he  Fourier  transformable  with  the  transform 

(..li  t  fl  I  COS  , 

II  -TTTTT- - TV-  (2  "> 

./(  I-P  I  A*(|  cos  c/>o ) 

Note  that  the  PO  basis  function  depends  on  the  incident  angle  Therefore, 


each  time  the  incident  angle  changes,  as  in  a  backscatter  pattern  computation, 
the  impedance  mat  rix  changes.  More  specific  ally,  as  0, (  changes,  the  rows  and 
columns  of  the  impedance  matrix  associated  with  the  1*0  basis  function  must  be 
recomputed. 

The  surface  wave  basis  function  is  also  used  when  Z$(x )  /  Z  over  a  semi¬ 
infinite  range.  However,  it  is  used  only  when  a  surface  wave  can  be  supported  by 
the  multilayered  slab  geometry.  Recall  that  if  a  surface  wave  exists  and  it  decays 
quickly  enough,  then  a  surface  wave  basis  function  does  not  need  to  be  included 
for  that  surface  wave.  The  surface  wave  basis  function  is  used  to  model  the  current 
variation  associated  with  a  surface  wave',  “launched”  from  the'  deviation  of  Zy(.r) 
from  Zc; Q,  which  travels  in  the  \  x  -  direction.  In  general,  more  than  one  surface 
wave  can  exist.  If  a  surface  wave  is  launched,  it  will  have  the  following  x  variation: 


whe, 


Refyy)  >  0 
1 111(75)  <  0 


(2.15) 


where  j 75  is  the  surface  wave  propagation  constant.  T  he  conditions  on  75  insure 
that,  the  surface  wave  decavs  as  it  travels  in  the  f  x  -  direction.  T  he  surface  wave 


basis  function  is  chosen  to  be 


n  if, 


:r  >  x  c 


>1  herwise 


(2.  Ifi) 


where  x$  is  typically  0  or  the  value  of  x  where  the  pulse  basis  functions  stop.  ,/, 

is  Fourier  t ransformable  with  the  transform 

,  7t  ■  t.s  >'\s- 

’/"  T0T~vi)' 


(2.17) 


To  use  surface  wave  basis  functions,  the  surface  wave  propagation  constants 
must  first  be  obtained.  However,  it  should  be  noted  that  surface  waves  will  not 
always  exist.  If  they  do  exist,  the  surface  wave  propagation  constants  for  a  mul¬ 
tilayered  slab  geometry  are  the  poles  of  the  function  F((3)  where  Z5 (|  is  replaced 
by  the  sheet  impedance  the  surface  waves  exist  on.  F(f3)  is  given  in  Appendix  A 
for  several  multilayered  slab  geometries.  In  general,  solving  for  the  surface  wave 
poles  will  result  in  a  transcendental  equation  which  can  be  solved  numerically.  The 
surface  wave  poles  must  satisfy  the  conditions  stated  in  Equation  (2.45).  When 
choosing  surface  wave  basis  functions,  those  associated  with  propagating  and  slowly 
decaying  surface  waves  should  be  chosen  first  since  they  are  the  dominant  modes. 

In  this  MM  solution,  surface  wave  basis  functions  are  included  only  in  the 
case  where  Zg{x)  is  entirely  in  free  space,  i.e.,  no  material  layers  exist.  However, 
surface  wfaves  cannot  exist  on  a  purely  real  sheet  impedance  in  free  space  [19].  In 
Figures  6(b)  and  6‘(c),  if  lm(Z.\)  ^  0,  then  a  surface  wave  will  travel  in  the  4  x  - 
direction  along  Zj\.  In  this  case,  solving  for  the  poles  of  Equation  (A. 20)  where 
Zsn  ZA>  the  surface  wave  propagation  constant  is  found  to  be 

g 2 

j IS  ~  F.jk[ ).  1  -  -A-  (2.48) 

N 

where  Zq  =  y/^1  is  the  characteristic  impedance  of  free  space.  The  root  of  Equa¬ 
tion  (2.48)  is  chosen  so  that  the  conditions  of  Equation  (2.45)  are  satisfied. 
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CHAPTER  III 


NUMERICAL  RESULTS 


This  chapter  presents  results  obtained  from  the  MM  procedure  outlined  in 
Chapter  II.  The  results  include  computed  backscatter  from  various  sheet  imped¬ 
ances  and  multilayered  geometries.  The  results  are  compared  with  measurements 
and  previously  calculated  results. 

3.1  Convergence  and  Current  Distribution  Results 

This  section  shows  convergence  data  on  the  echo  width  scattered  by  a  perfectly 
conducting  half  plane  for  edge-on  incidence.  Figure  7  shows  the  edge-on  echo  width 
versus  the  width  of  the  pulse  basis  function  region,  labeled  W  in  the  figure.  The 
PO  basis  function  is  required.  The  frequency  is  300  MHz  so  that  A  =  1  meter. 
Each  pulse  basis  function  is  0.1A  wide  and  the  W  is  varied  from  0  to  3A.  Tims,  the 
number  of  pulses  varies  from  0  to  30  pulses.  It  can  be  seen  that  the  echo  width 
converges  to  within  0.25  dB  of  the  exact  value  at  a  pulse  region  width  of  about 
1A.  This  appears  to  be  the  best  this  MM  solution  can  do  for  a  pulse  basis  function 
width  of  0.1A.  If  greater  accuracy  is  desired,  a  smaller  pulse  width  can  be  used. 

Figure  8  shows  the  edge-on  echo  width  versus  the  width  of  the  individual  pulse 
basis  functions,  labeled  I)  in  the  figure.  Once  again,  the  l’()  basis  fund  ion  is  used 
and  the  frequency  is  300  MHz.  The  pulses  are  placed  over  the  first  wavelength  (YV 
=  A)  and  D  is  varied  from  0.033  to  0.25.  At  D  — 0.033  there  are  30  pulses  and  at 


28 


D=0.25  there  are  4  pulses.  The  MM  solution  converges  nicely,  thus  demonstrating 
the  accuracy  of  the  solution. 

The  MM  current  distribution  for  a  typical  basis  function  expansion  is  com¬ 
pared  to  the  exact  current  distribution  on  the  PEC  half  plane  [5]  in  Figure  9.  The 
MM  current  distribution  uses  a  PO  basis  function  starting  at  xpQ  --  0.  Also, 
10  pulse  basis  functions  of  width  0.1  A  each  are  used  over  the  first  wavelength  of 
the  PEC  half  plane.  The  frequency  is  300  MHz  (A  =  1  meter)  and  the  current 
distribution  is  plotted  over  the  first  two  wavelengths  of  the  half  plane.  The  MM 
and  exact  currents  are  in  reasonable  agreement  with  one  another.  For  the  range 
0  <  x  <  A,  the  current  consists  of  pulses  superimposed  on  the  PO  current,  thus 
giving  the  distribution  the  slightly  modified  staircase  shape.  For  x  >  A  the  current 
distribution  consists  only  of  the  PO  current.  This  can  be  seen  in  the  constant 
magnitude  and  linear  phase  of  the  current  in  this  region.  It  is  interesting  to  note 
that  the  MM  current  at  the  center  of  the  first  pulse  is  about  1.44  times  the  exact 
current  at  this  location.  The  exact  current  has  a  dependance  near  the  edge  of 
the  half  plane.  The  average  value  of  the  exact  current  over  the  first  pulse  width 
divided  by  its  value  at  the  center  of  the  first  pulse  is  y/2  ~  1.41.  Thus,  despite  the 
appearance  of  the  current  plot,  the  MM  current  distribution  does  a  very  good  job 
of  approximating  the  exact  current  over  the  first  pulse  region. 

3.2  Backscatter  Versus  Frequency  Results 

This  section  presents  a  comparison  of  measured  backscatter  from  a  sheet  im¬ 
pedance  located  in  a  A  —  2  layered  slab  with  results  obtained  from  this  MM 
solution.  Measured  and  computed  echo  width  versus  frequency  are  compared. 
The  MM  computations  were  obtained  using  an  impedance  matrix  interpolation 
method  [20].  The  inserts  in  Figures  10  and  11  show  the  multilayered  slab  geome- 
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ECHO  WID' 


W/A  =  width  of  pulse  region 


Figure  7:  Convergence  curve  showing  edge-on  echo  width  versus  pulse  region 

width  for  P FC  half  plane. 
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D/A=  width  of  each  pulse 


Figure  8:  Convergence  curve  showing  edge-on  echo  width  versus  individual  pulse 

width  for  PEC  half  plane. 
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Figure  9:  Comparison  of  MM  current  with  exact,  current  for  edge-on  incidence  to 

a  P EC  half  plane. 
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try  analyzed  here. 

The  geometry  consists  of  a  sheet  impedance  like  that  shown  in  Figure  3(b) 
given  by  Equation  (2.2b  Referring  to  Equation  (2.2)  the  sheet  impedance  is  defined 
by  —  650  and  Z  \  0  or  oo.  Tlie  finite  width  where  Z  4  exists  is  L  ~  7.62cm. 

This  sheet  impedance  is  located  in  a  K  =  2  layered  slab.  Both  layers  are  lossless 
and  have  permeability  //q  arid  relative  dielectric  constant  er  —  2.6.  The  thickness 
of  the  layers  is  T\  -  75  -  0.57cm. 

Figures  10  and  11  show  plots  of  measured  results  compared  with  results  ob¬ 
tained  from  this  MM  solution.  Figure  10  corresponds  to  Z\  —  oo  and  Figure  11 
corresponds  to  Z ^  —  0.  The  figures  show  backscatter  versus  frequency  at  </> q  —  30°. 
The  MM  computations  were  made  using  impedance  matrix  interpolation.  The  im¬ 
pedance  matrix  was  computed  every  A/  =  1  GHz  and  echo  width  was  calculated 
every  A /  —  25  MHz.  33  pulse  basis  functions  were  used  in  the  range  from  4  to  8 
GHz  and  50  pulse  basis  functions  were  used  from  8  to  12  GHz.  This  corresponds 
to  a  pulse  width  of  about  O.lAf)  at  the  highest  frequency  in  each  range  where  \p 
is  the  wavelength  in  the  dielectric  layer.  The  solid  curves  are  the  calculated  data 
and  the  dashed  curves  are  the  measured  data.  The  results  agree  reasonably  well, 
especially  for  Z -  00.  The  data  took  about  22  minutes  per  plot  to  compute  on 
a  VAX  8550. 

3.3  Backscatter  From  a  Sheet  Impedance  in  Free  Space 

This  section  presents  backscatter  from  a  sheet  impedance  half  plane  in  free 
space.  The  geometry  is  sketched  in  the  inserts  of  Figures  12  and  13.  Several 
values  for  purely  resistive  and  purely  reactive  sheet  impedances  are  considered 
with  backscatter  patterns  of  their  echo  width  included. 

Figure  12  shows  the  computed  backscatter  echo  width  of  a  purely  resistive 
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Figure  10:  Backscatter  by  a  shoe!  impedance  located  in  a  two  layered  slab. 

Z511  — -  0511  and  7j  \  ~  00. 
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sheet  impedance  half  plane  in  free  space.  The  sheet  impedance  is  of  the  form 


shown  in  Figure  3(c)  given  in  Equation  (2.3)  where  Zg u  oo  and  Z.j  -  Rg 

with  Rg  —  0,100,500.2000,10*12.  The  echo  width  is  given  in  dB-m  and  was 
calculated  at  300.0  MHz.  There  were  10  pulse  basis  functions  of  width  0.1  A 
10cm  each  covering  the  first  wavelength  of  the  half  plane,  furthermore,  the  PO 
basis  function  (starting  at  r pQ  -  0)  had  to  be  included,  making  a  total  of  II 
basis  functions.  No  surface  wave  basis  function  was  used  since  a  purely  real  sheet 
impedance  cannot  support  a  surface  wave.  The  edge-on  backseat  ter  at  180° 

is  in  excellent  agreement  with  values  previously  calculated  by  Senior  [5j.  These 
points  are  indicated  as  large  dots  at  180°  in  figure  12.  Also,  the  MM  results 
computed  here  check  out  very  well  against  the  asymptotic  results  of  Ly  [15].  Plots 
comparing  these  results  were  too  similar  to  include  here.  The  data  in  Figure  12 
took  about  28  minutes  of  CPU  time  to  compute  on  a  VAX  8550. 

Similarly,  Figure  13  shows  the  computed  backscatter  echo  width  of  a  purely 
reactive  sheet  impedance  half  plane  in  free  space.  The  sheet  impedance  is  of 
the  same  form  given  in  Equation  (2.3)  where  ZgQ  -  oo  and  Z -j Xg  with 
A* 5  =  0,100,500,2000,  10'*fi.  The  echo  width  is  given  in  dB-m  and  was  calculated 
at  300.0  MIIz.  There  were  10  pulse  basis  functions  of  width  0.1A  =  10cm  over  the 
first  wavelength  of  the  half  plane.  I'he  imaginary  sheet  impedance  can  support 
a  surface  wrave  so  the  surface  wave  basis  function  was  included.  It  started  at 
xg  =  0  and  had  propagation  constant  given  by  Equation  (2.48).  With  the  PO  basis 
function,  a  total  of  12  basis  functions  were  used.  The  imaginary  sheet  impedance 
scattering  results  also  check  out  very  well  with  the  work  of  Ly  [15].  The  data  in 


Figure  13  took  about  31  minutes  to  compute  on  a  VAX  8550. 
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Figure  12:  Backseat  ter  by  a  purely  real  sheet  impedance  half  plane  in  free  space 
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Figure  13:  Backscatter  by  a  purely  imaginary  sheet  impedance  half  plane  in  free 

space. 


3.4  Backscatter  From  a  Tapered  Sheet  Impedance  Half  Plane  in  Free 
Space 

This  section  presents  results  of  backscatter  from  a  tapered  sheet  impedance 
half  plane  in  free  space.  The  geometry  is  sketched  in  the  insert  of  Figure  14  and 
the  sheet  impedance  variation  is  described  below.  Backscatter  patterns  of  echo 
width  are  included  for  several  different  linear  impedance  tapers.  Also,  backscatter 
patterns  from  a  linearly  tapered  and  a  exponentially  tapered  sheet  impedance  are 
compared  against  results  from  Newman  [8]  [9]. 

The  sheet  impedance  variation  is  similar  to  that  shown  in  Figure  3(d)  given 
by  Equation  (2.4).  The  sheet  impedance  tapers  linearly  from  lOOOfi  to  014  over 
the  width  L.  For  x  >  L  the  sheet  impedance  is  zero,  i.e.,  a  perfectly  conducting 
half  plane.  Referring  to  Equation  (2.4),  the  sheet  impedance  considered  here  is 
given  as 

oo  if  x  <  0 

Zs(x)  =  '  1000  (l-f)  if  0  <  x  <  L  (3-1) 

0  if  x  >  L. 

Backscatter  patterns  were  computed  for  L  —  0,  A/4,  A/2,  A  and  2A.  The  backscatter 
patterns  of  echo  width  in  dB-m  are  shown  in  Figure  14  and  were  computed  at  300 
MHz.  In  each  case,  pulse  basis  functions  of  width  0.05A  were  placed  over  the 
impedance  taper  region  and  the  first  wavelength  of  the  PEC,  i.e.,  they  were  used 
from  x  =  0tox  =  Z/  +  A.  The  PO  basis  function  had  to  be  included  and  it  started 
at  xpQ  =  0.  The  data  of  Figure  14  took  about  145  minutes  of  CPU  time  on  a 
VAX  8550. 

For  comparison,  backscatter  patterns  from  linearly  tapered  and  exponentially 
tapered  sheet  impedances  were  compared  with  Newman’s  results.  The  width  of 
the  impedance  taper  region  is  L  —  A  in  both  cases.  The  linear  impedance  taper 
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is  that  of  Equation  (3.1).  The  exponential  impedance  tapers  from  1000Q  to  1 00 
exponentially  and  is  given  as 

oc  if  x  <  0 

Zs(x)  =  i0e-4.61(*-L)/A  if  o  <x<L  (3.2) 

0  if  x  >  L. 

Figure  15  shows  the  results  of  this  MM  solution  compared  with  Newman’s  results 
for  both  impedance  tapers.  The  results  are  in  good  agreement,  thus  demonstrating 
the  accuracy  of  this  solution. 

3.5  Backscatter  From  a  Tapered  Sheet  Impedance  Discontinuity  in 
Free  Space 

This  section  presents  results  of  backscatter  from  a  tapered  sheet  impedance 
discontinuity  in  free  space.  The  geometry  is  sketched  in  the  insert  of  Figure  16. 
The  problem  considered  in  this  section  is  similar  to  that  of  the  previous  section 
with  the  main  difference  being  that  the  “background”  impedance  Zgo  is  100Q 
instead  of  infinity.  Backscatter  patterns  are  included  for  several  linear  impedance 
tapers.  Also,  backscatter  patterns  from  a  linearly  tapered  and  a  cosinusoidally 
tapered  sheet  impedance  are  compared. 

The  sheet  impedance  tapers  linearly  from  lOOfi  to  Off  over  the  width  L  and 
can  be  given  by 

100  if  x  <  0 

ZS{x)  =  ■  100  (l  -  f )  if  0  <  x  <  L  (3.3) 

0  if  x  >  L. 

Backscatter  patterns  were  computed  for  L  =  0,  A/4,  A/2,  A  and  2A  at  300  MHz  and 
are  shown  in  Figure  16.  In  each  case,  pulse  basis  functions  of  width  0.05A  were 
placed  from  x  =  0to*  —  L  +  A.  The  PO  basis  function  had  to  be  included  starting 
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Figure  15:  Comparison  of  backscatter  by  linearly  and  exponentially  tapered  sheet 

impedance  half  planes  in  free  space. 
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at  xpo  —  0.  The  data  of  Figure  14  took  about  145  minutes  of  CPU  time  on  a 
VAX  8550. 

A  cosinusoidal  sheet  impedance  taper  from  1000  to  Of 1  of  width  L  =  A  was 
also  analyzed.  The  cosinusoidal  impedance  taper  is  given  as 


Zg(x)  =  50  +  50  cos 


for  0  <  x  <  L. 


(3.4) 


The  backscatter  from  this  impedance  taper  is  contrasted  with  the  linear  imped¬ 
ance  taper  of  the  same  width  in  Figure  17.  Both  sheet  impedances  had  the  same 
basis  function  expansion  of  40  pulse  basis  functions  of  width  0.05A  over  the  first  2 
wavelengths.  Both  cases  also  used  the  PO  basis  function. 
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Figure  16:  Backscatter  by  a  linearly  tapered  sheet  impedance  discontinuity  in 

free  space. 
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Figure  17:  Comparison  of  backscatter  by  linearly  and  cosinusoidally  tapered 
sheet  impedance  discontinuities  in  free  space. 
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CHAPTER  IV 


DESCRIPTION  OF  COMPUTER  CODE 

This  chapter  describes  briefly  a  computer  code  which  implements  the  MM 
solution  outlined  in  Chapter  II.  The  computer  code  has  been  named  the  uZs(x)” 
Code  (ZSX)  by  the  author.  ZSX  can  analyze  the  three  multilayered  slab  geome¬ 
tries  considered  in  Appendices  A  and  B.  These  geometries  are  a  sheet  impedance 
located  in  a  K  =  2  layered  slab,  on  a  material  half  space,  and  in  free  space.  ZSX 
can  calculate  bistatic  or  backscatter  patterns,  current  distributions,  and  frequency 
sweep  data.  The  inputs  and  outputs  for  ZSX  are  explained  in  this  chapter. 

4.1  Input  and  Output  Files 

ZSX  utilizes  one  input  file  and  two  output  files.  The  input  file,  INF1L.DAT, 
contains  all  the  input  data  to  ZSX,  i.e.,  all  the  input  data  are  read  from  IN- 
FIL.DAT.  INFIL.DAT  must  be  assigned  to  logical  unit  8.  The  first  output  file, 
OUTFL.DAT,  contains  the  output  of  ZSX.  After  a  run,  OUTFL.DAT  contains 
the  input  data  used  in  that  run,  the  problem  geometry  and  basis  function  expan¬ 
sion  used,  and  a  tabular  listing  of  the  electromagnetic  calculations  if  any  were 
made.  OUTFL.DAT  must  be  assigned  to  logical  unit  9.  The  second  output  file. 
PTPLOT.DAT,  contains  tabular  listings  of  the  electromagnetic  calculations.  PT- 
PLOT.DAT  is  intended  for  plotting  purposes.  PTPLOT.DAT  will  contain  either 
a  bistatic  or  backscatter  pattern,  a  current  distribution,  or  frequency  sweep  data. 
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c 

C  READ  1 

READ (8, *) NGO, IWRZM, IWRC, IGMTYP , IMF, IMS, IMP,NSH, BTMAX, ICALC 
C 

C  READ  2 

READ ( 8 , * ) F 
C 

C  READ  3 A 

IF ( IGMTYP .EQ.2) READ (8 ,  *) ER1 , TDE1 , UR1, TDM1 
C 

C  READS  3B  AND  3C 

IF ( IGMTYP . EQ . 3 ) THEN 

READ (8, *) ER1 , TDE1 , UR1 , TDM1 , T1 
READ ( 8 , * ) ER2 , TDE2 , UR2 , TDM2 , T2 
ENDIF 
C 

C  READ  4A 

READ (8, *) ZSO, ZS 
C 

C  READ  4B 

READ {8, *) ITP, ZA, ZB, WD1 
C 

C  READ  5 

READ < 8 , * ) WD2 , SGXW , XI , X2 
C 

C  READ  6A 

IF ( ICALC . EQ . 0 ) READ ( 8 , * ) IPAT , DPHI , PHBST 
C 

C  READ  6B 

IF ( (ICALC. EQ.l) .OR. (ICALC. EQ.2) ) READ ( 8 , * ) FMC1 , FMC2 , DFZ, DFC, 
&PHIN , PHSC 

C 

C  READ  6C 

IF (IPAT. EQ. 3) READ ( 8 , * ) XI , XF , XS 


Figure  18:  The  FORTRAN  READ  statements  in  ZSX. 
PTPLOT.DAT  must  be  assigned  to  logical  unit  10. 

4.2  Input  Data 

The  input  data  to  ZSX  are  explained  in  this  section.  The  input  data  are 
used  to  describe  to  the  program  the  problem  geometry  and  indicate  the  desired 
electromagnetic  calculations.  The  input  data  are  obtained  via  FORTRAN  READ 
statements.  The  READ  statements  and  the  parameters  defined  by  them  will  be 
explained  in  this  section.  Figure  18  shows  all  the  READ  statements  contained  in 
ZSX.  Note  that  not  every  READ  statement  will  be  executed  in  every  program  run. 
The  input  file,  INFIL.DAT,  contains  the  data  to  be  read  by  the  READ  statements. 
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4.2.1  READ  1:  Run  Control  Parameters 

READ  1  defines  the  following  run  control  parameters: 

NGO  -  run  indicator. 

—  0  implies  input  and  print  out  problem  description  and  then  stop,  i.e.,  do 
not  make  any  electromagnetic  calculations. 

=  1  implies  input  problem  description  and  then  perform  the  desired  electro¬ 
magnetic  calculations.  An  NGO  =  0  run  should  precede  an  NGO  -=  1  run  as 
this  allows  the  user  to  verify  the  accuracy  of  the  problem  description  defined 
in  the  input  file,  INFIL.DAT. 

IWRZM  =  indicator  for  writing  the  impedance  matrix  to  the  output  file. 

—  0  implies  do  not  write  the  impedance  matrix  to  the  output  file. 

■=  1  implies  write  the  impedance  matrix  to  the  output  file.  Note  that  this 
will  result  in  N2  lines  of  output.  Furthermore,  the  impedance  matrix  will  be 
printed  at  every  angle  if  a  physical  optics  (PO)  basis  function  is  used  and  a 
backscatter  pattern  is  defined. 

IWRC  =  indicator  for  writing  out  both  the  voltage  and  current  vector. 

—  0  implies  do  not  write  out  the  voltage  or  current  vector. 

1  implies  write  out  both  the  voltage  and  current  vector.  Note  that  for 
backscatter  patterns  this  will  result  in  the  voltage  and  current  vector  being 
printed  at  every  angle. 

IGMTYP  =  indicator  for  the  type  of  multilayered  slab  geometry. 

—  1  implies  a  sheet  impedance  in  free  space. 

—  2  implies  a  sheet  impedance  on  a  material  half  space. 

=  3  implies  a  sheet  impedance  located  in  a  K  -  2  layered  slab. 


■18 


IMF  =  indicator  for  using  a  physical  optics  (PO)  basis  function. 

=  0  implies  do  not  use  a  PO  basis  fu  on.  Use  this  in  the  case  that 
Zg(x)  Zg o  over  a  finite  range  of  x. 

=  1  implies  use  a  PO  basis  function.  Use  this  in  the  case  that  Zg(x)  ^  Zg q 
over  a  semi-infinite  range  of  x. 

IMS  =  indicator  for  using  a  surface  wave  basis  function. 

=  0  implies  do  not  use  a  surface  wave  basis  function. 

=  1  implies  use  a  surface  wave  basis  function.  IMS  should  be  set  equal  to 
1  only  if:  1)  the  sheet  impedance  is  NOT  purely  real  and  is  in  free  space 
(IGMTYP  =  1),  and  2)  the  range  where  Zg(x)  ^  ZgQ  is  semi-  infinite.  ZSX 
allows  for  only  one  surface  wave  basis  function  for  the  case  that  the  sheet 
impedance  is  in  free  space.  ZSX  computes  the  surface  wave  propagation 
constant  according  to  Equation  (2.48). 

IMP  =  indicator  for  using  pulse  basis  functions. 

=  0  implies  do  not  use  any  pulse  basis  functions. 

=  1  implies  use  pulse  basis  functions.  This  will  almost  always  be  the  case. 

NSH  =  number  of  Simpson’s  rule  integration  segments  used  per  half  cycle  of 
the  oscillatory  part  of  the  integrand  of  the  spectral  integral  expression  for 
Zmn  (Equation  (2.31).)  NSH  =  4  typically  but  can  be  increased  for  greater 
accuracy. 

BTMAX  =  upper  and  lower  limits  of  integration  used  in  the  expression  for  Zmn 
in  terms  of  the  maximum  wavenumber  of  the  multilayered  slab  materials, 
i.e.,  limits  of  integration  for  Equation  (2.31)  are  iBTMAXxfcp  where 
k[)max  is  the  maximum  wavenumber  of  any  material  in  the  multilayered  slab. 
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BTMAX  must  be  great  enough  to  allow  for  convergence  of  the  integral.  H  I  - 
MAX  =  50.0  typically.  To  check  for  convergence,  BTMAX  can  be  increased 
until  the  impedance  matrix  elements  are  stable. 

ICALO  -  indicator  for  performing  frequency  sweep  computations. 

-=  0  implies  do  not  perform  a  frequency  sweep  computation. 

—  1  implies  perform  a  frequency  sweep  computation  using  standard  quad¬ 
ratic  interpolation. 

=  2  implies  perform  a  frequency  sweep  computation  using  improved  quad¬ 
ratic  interpolation.  This  should  never  be  done  if  a  HO  basis  function  is  used 
(IMF  =  1)  or  if  a  surface  wave  basis  function  is  used  (IMS  1). 

4.2.2  READ  2:  Frequency 

READ  2  defines  F  the  frequency  in  MHz  at  which  any  pattern  or  current 
distribution  calculations  are  made.  If  a  frequency  sweep  computation  is  desired, 
then  the  frequencies  are  determined  via  READ  6B.  However,  READ  2  will  be 
executed  in  either  case. 

4.2.3  READ  3:  Material  Parameters 

Note  that  READ  3A  will  be  executed  only  if  IGMTY  1*  2  and  that  READ 

3B  and  3C  will  be  executed  only  if  IGMTY  P  —  3  (see  READ  I).  READ  3A 
determines  the  material  parameters  of  the  half  space  if  IGMTYP  -  2. 

ER]  —  relative  real  part  of  permittivity  in  material  half  space. 

TDE1  —  loss  tangent  of  permittivity  in  material  half  sparr. 

UR!  =  relative  real  part  of  permeability  in  material  half  space. 


TDM1  -  loss  tangent  of  permeability  in  material  half  space. 


HEAD  3B  and  3C  determine  the  material  parameters  of  the  two-layered  slab 
if  IGMTYP  =  3. 

ER1  =  relative  real  part  of  permittivity  in  layer  k  —  1. 

TDEl  =  loss  tangent  of  permittivity  in  layer  k  —  1. 

URl  =  relative  real  part  of  permeability  in  layer  k  —  1. 

TDM1  =  loss  tangent  of  permeability  in  layer  k  —  1. 

Tl  =  thickness  in  meters  of  layer  k  —  1. 

ER2  =  relative  real  part  of  permittivity  in  layer  k  =  2. 

TDE2  =  loss  tangent  of  permittivity  in  layer  k  =  2. 

UR2  =  relative  real  part  of  permeability  in  layer  k  =  2. 

TDM2  =  loss  tangent  of  permeability  in  layer  k  —  2. 

T2  =  thickness  in  meters  of  layer  k  =  2. 

4.2.4  READ  4:  Sheet  Impedance 

READ  4A  defines  the  constant  values  of  the  sheet  impedance  Zg(x)  as  follows: 

ZSO  =  “background”  constant  sheet  impedance  Zg o-  If  the  user  is  interested  in 
a  case  where  ZgQ  — »  oo,  such  as  scattering  from  a  resistive  half  plane  in  free 
space,  then  simply  set  ZSO  equal  to  a  very  large  number.  Experience  has 
shown  that  ZSO  =  10^  gives  good  results  in  these  cases. 

ZS  =  constant  value  of  the  sheet  impedance  (not  Zg q)  in  the  range  where  the 
sheet  impedance  is  constant,  but  not  equal  to  Zg q.  For  example,  set  ZS  = 
ZA  for  the  sheet  impedances  shown  in  Figures  3(b),  3(c)  and  3(d). 
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The  width  where  Zg(x)  is  constant  (not  Zg q)  is  taken  to  he  finite  if  only  pulse 
basis  functions  are  used  (see  Figure  6(a).)  This  width  is  taken  to  be  semi-infinite 
going  to  infinity  in  the  +£■  -  direction  if  a  PO  or  surface  wave  basis  function  is 
used  (see  Figures  6(b)  and  6(c).) 

READ  4B  can  be  used  to  define  an  impedance  taper  region.  An  impedance 
taper  can  account  for  a  non-constant  sheet  impedance  like  that  shown  in  Figure 
3(d).  ZSX  can  allow  for  only  one  sheet  impedance  taper.  READ  4B  defines  the 
following: 

ITP  =  indicator  for  an  impedance  taper. 

=  0  implies  do  not  include  an  impedance  taper. 

=  1  implies  include  a  constant  impedance  taper. 

=  2  implies  include  a  linear  impedance  taper. 

=  3  implies  include  an  exponential  impedance  taper. 

=  4  implies  include  a  cosine  impedance  taper. 

ZA  =  impedance  value  at  the  start  of  the  impedance  taper.  ZA  will  be  the  value 
of  a  constant  impedance  taper  if  ITP  =  1.  Note  that  ZA  input  via  this 
READ  statement  is  not  to  be  confused  wHh  Z of  Figures  3  and  6. 

ZB  —  impedance  value  at  the  end  of  the  impedance  taper. 

WD1  =  width  in  meters  of  the  impedance  taper  region. 

If  ITP  >  0  then  there  will  be  an  impedance  taper  region.  This  impedance  taper 
region  starts  at  x  —  — WD1  and  extends  to  ,r  0.  If  ITP  J  then  (he  sheet 
impedance  will  have  a  constant  value  of  ZA  over  (his  entire  range  of  x.  If  ITP  - 
2,3  or  4  then  the  sheet  impedance  will  have  value  ZA  at  x  -WD1  and  value 
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ZB  at  x  =  0.  In  these  cases,  the  sheet  impedance  will  taper  from  ZA  to  ZB  either 
linearly,  exponentially  or  cosinusoidally,  depending  on  ITP. 

4.2.5  READ  5:  Basis  Function  Description 

READ  5  defines  the  layout  of  the  basis  function  expansion  through  the  use  of 
the  parameters  explained  below. 

WD2  =  width  of  constant  ZS  where  pulse  basis  functions  will  be  used.  This 
parameter  has  meaning  only  if  IMP  =  1.  If  IMP  =  1  then  pulse  basis 
functions  are  placed  from  x  —  0  to  x  =  WD2. 

SGXW  =  maximum  segment  size  of  pulse  basis  functions  in  minimum  wavelengths 
of  a  layer.  SGXW  should  not  exceed  0.25  but  has  typically  been  chosen  as 
0.1  or  less. 

XI  =  the  x  value  in  meters  where  the  physical  optics  basis  function  starts.  XI 
has  meaning  only  if  IMF  =  1.  XI  is  analogous  to  xpQ  of  Equation  (2.13) 
and  is  typically  chosen  as  0.0  or  WD2. 

X2  =  the  x  value  in  meters  where  the  surface  wave  basis  function  starts.  X2  has 
meaning  only  if  IMS  -  1.  X2  is  analogous  to  xg  of  Equation  (2.16)  and  is 
typically  chosen  as  0.0  or  WD2. 

Note  that  if  an  impedance  taper  region  is  specified  (ITP  >  0  in  READ  4B)  then 
pulse  basis  functions  are  automatically  placed  over  this  region,  i.e.,  if  ITP  >  0 
then  pulse  basis  functions  are  placed  from  x  --  —  VYD1  to  x  =  0.  For  pulse  basis 
functions  to  be  placed  over  the  impedance  taper  region,  IMP  need  not  be  set  be  set 
to  1  in  READ  1.  The  width  of  these  pulse  basis  functions  over  WDl  is  determined 
by  SGXW. 
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4.2.6  READ  6:  Electromagnetic  Calculations 

Note  that  READ  6A  is  executed  only  if  ICAEC  --  0  and  READ  611  is  executed 
only  if  ICAEC  —  1  or  2  (see  READ  1).  Similarly,  REA  D  6C  is  executed  only  if  1  PAT 
—  3  in  READ  6A.  READ  6A  specifies  whether  a  pattern  or  a  current  distribution  is 
to  be  computed.  It  also  defines  the  pattern  using  the  parameters  explained  below. 

IPAT  =  indicator  for  computing  either  a  pattern  or  a  current  distribution. 

—  0  implies  do  NOT  compute  a  pattern  or  a  current  distribution. 

=■  1  implies  compute  a  backseat. ter  pattern. 

=  2  implies  compute  a  bistatic  scattering  pattern. 

■=  3  implies  compute  a.  current  distribution  (see  READ  6C). 

DPHI  =  the  angle  increment,  for  pattern  computations.  If  IPAT  —  1  or  2  the 
scattering  pattern  will  be  computed  every  DPI1I  degrees. 

PHBST  =  the  angle  of  the  incident  plane  wave  for  bistatic  scattering  patterns  or 
current  distributions. 

READ  6B  defines  a  frequency  sweep  computation  using  the  parameters  ex¬ 
plained  below. 

FMCl  —  beginning  frequency  in  MIIz  for  a  frequency  sweep  computation. 

EMG'2  —  ending  frequency  in  MHz  for  a  frequency  sweep  computation. 

DEZ  —  frequency  step  size  in  MHz  for  calculating  the  impedance  matrix. 

DEC  —  frequency  step  size  in  MHz  for  calculating  the  scattered  field. 

PIIIN  angle  of  the  incident  field  in  degrees. 

PHSC  angle  of  the  scattered  field  in  degrees. 
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If  standard  interpolation  is  used  (ICALC  =  1  in  READ  1)  then  the  frequency  step 
size  DFZ  in  MHz  is  typically  set  to  A/ji// 2  where 


A/a/  = 


300 

2 T 


MHz. 


Here  L  is  the  width  over  which  pulse  basis  functions  are  placed.  Assuming  that 
pulse  basis  functions  are  used  (IMP  =  1  in  READ  1),  then  L  —  VVD2  if  no  imped¬ 
ance  taper  is  used  (ITP  =  0  in  READ  4B)  and  L  =  WD1  4-  WD2  if  an  impedance 
taper  is  used  (ITP  >  0  in  READ  4B).  If  improved  interpolation  is  used  (ICALC  — 
2  in  READ  1)  then  DFZ  can  be  increased  to  or  more.  Recall  that  improved 

interpolation  should  not  be  used  if  either  a  PO  or  surface  wave  basis  function  is 
used. 

READ  6C  defines  the  current  distribution  via  the  parameters  explained  below. 


XI  =  the  initial  value  of  x  in  meters  for  the  current  distribution. 


XF  =  the  final  value  of  x  in  meters  for  the  current  distribution. 


XS  =  the  step  size  of  x  in  meters  for  the  current  distribution. 

If  a  current  distribution  is  specified  (IPAT  =  3  in  READ  6A)  then  the  current 
induced  by  a  plain  nave  incident  from  ~  PHBST  is  computed.  The  current 
distribution  is  computed  from  x  =  XI  to  x  =  XF  at  increments  of  XS. 
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CHAPTER  V 


SUMMARY 


This  report  has  described  the  integral  equation  and  MM  solution  for  two- 
dimensional  TM  scattering  by  a  variable  sheet  impedance  in  a  multilayered  slab. 
An  integral  equation  was  derived  for  the  unknown  surface  current  on  the  scat 
tering  portion  of  the  sheet  impedance.  The  moment  method  was  applied  to  the 
integral  equation  to  obtain  an  approximate  expansion  for  the  unknown  surface 
current.  Simple  expressions  for  the  impedance  matrix  and  voltage  vector  elements 
used  in  the  MM  solution  were  obtained.  The  far  field  of  the  basis  functions  were 
determined  asymptotically  so  the  far  zone  scattered  field  could  be  obtained  analyt¬ 
ically.  Numerical  results  were  presented  and  some  results  were  compared  against 
measured  or  previously  calculated  results. 

One  important  feature  of  this  MM  solution  is  that  the  impedance  matrix 
elements  and  the  far  zone  fields  depend  upon  spectral  functions  specific  to  the 
geometry  of  the  multilayered  slab.  Thus,  different  spectral  functions  can  be  used 
to  apply  the  MM  solution  to  different  multilayered  slab  geometries.  Spectral  func¬ 
tions  were  obtained  for  a  sheet  impedance  located  in  a  A'  2  layered  slab,  on 
a  material  half  space,  or  in  free  space.  A  user-oriented  computer  code  was  writ¬ 
ten  to  implement  this  MM  solution  and  can  account  for  the  three  geometries  just 
mentioned.  A  user’s  manual  for  this  code  was  given  in  Chapter  IV. 

Future  work  in  this  area  might  include  the  computation  of  transverse  electric 
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(TE)  scattering  by  the  same  geometries  analyzed  here.  This  solution  will  result 
in  different  spectral  functions  and  will  require  different  basis  function  expansions. 
Furthermore,  the  volume  equivalence  theorem  applied  to  the  thin  dielectric  slab 
may  not  be  valid  for  near  edge-on  incidence,  and  this  problem  may  have  to  be 
treated  separately. 

Another  idea  for  further  study  is  a  more  complete  analysis  of  the  surface  wave 
poles  and  the  surface  wave  basis  functions.  The  MM  solution  presented  in  this 
report  only  used  surface  wave  basis  functions  when  the  sheet  impedance  was  in  free 
space.  If  the  dominant  surface  wave  poles  can  be  determined  for  more  complicated 
layered  geometries,  then  their  corresponding  surface  wave  basis  functions  can  be 
included  in  the  MM  current  expansions  for  these  geometries. 
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APPENDIX  A 


THE  DETERMINATION  OF  F(Q)  FOR  SEVERAL 
MULTILAYERED  SLAB  GEOMETRIES 

This  appendix  derives  the  function  F(/?)  of  Equation  (2.29)  for  the  two  layered 
geometry  shown  in  Figure  19.  Equations  (2.22)  and  (2.25)  are  used  to  express  the 
electric  field  and  the  tangential  magnetic  field.  The  conditions  on  these  fields, 
listed  in  Section  2.4,  are  enforced  and  the  function  F(/?)  is  solved  for.  From  this 

result,  the  function  F(fl)  is  also  found  for  the  two  special  limiting  cases  shown  in 

Figure  20. 

The  electric  field  is  expressed  as 

Bit*,*)  /IOO[  AU1)e  ***  |  »*(,»),’**  |  (.11) 

7-00  1  1 

and  the  tangential  magnetic  field  is  expressed  as 

«„(*, *)  =  —  /'”lt  r"lrdi)  (,1.2) 

JVftk  J  °° 

for  k =0, 1 ,2,3.  7 ^  is  given  by 

kk  ~T  ^x'f'kn 

Ik  H  ",,ere  '  Re(yA.)  >  0  (  l-d) 

•'»(■>*:)  °- 

The  first  condition  on  the  fields  is  the  radiation  condition  as  r  -■*  oo.  This 
implies  that  there  are  only  upward  travelling  waves  in  region  k  0  and  only 
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VA~  o  r  °r 

region:  k.  =  0 


Figure  19:  Sketch  of  two  layered  geometry  containing  sheet  impedance  Zs(x). 

downward  travelling  waves  in  region  k  =  3.  To  enforce  these  properties,  it  must 
be  true  that 


Bi(0)  =  0  and  Al((3)  =  0.  (AA) 

The  second  condition  is  continuity  of  the  tangential  electric  field  across  each 
interface.  Enforcing  this  condition  yields  the  following: 


Aim  e"7oTl 


Aln[0)e 


+Bln(P)e^T'}  =  0 


(AS) 

(A- 6) 

(.4.7) 

The  third  condition  is  continuity  of  the  tangential  magnetic  field  across  each 
interface  except  the  z  —  0  interface.  Enforcing  this  condition  yields  the  following: 


A\V>)  +  B^P)]  -  \Al(0)  +  Bl(0)\  =0 
A2n(0)e^  +  BldS)c-^  \  -  [B’08)e->»rs|  =  0. 


70 

JW/UQ 


A°n(f3)e~^Ti  +-2!-\Al(0)e-7iTl_B\(j3)e7iT1]=o  {A8) 

111)  II.  1 


59 


72 


A2n(0)eT*r2  -  B*(0)e  ^T2 
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Blm  e  73?2  -  0.  (,1.9) 


JUV2  1  1  JUfi  3 

The  fourth  and  final  condition  is  discontinuity  of  the  tangential  magnetic  field, 

bv  the  total  surface  current  produced  by  across  the  z  -  0  interface,  i.e., 


zx  ihS-  H^1 


J 


(•‘•10) 


where  k  and  k  -f  1  represent  the  regions  immediately  above  and  below  the  2  -  0 
interface  respectively,  and  is  the  total  surface  current  produced  by  Jn.  Since 
contains  only  a  y  -  directed  component,  then  for  the  K  —  2  layered  slab  this 
condition  becomes 

II L  -  4  =  Jn-  (-4.M) 

The  ( y  -  polarized)  surface  current  produced  by  Jn,  at  the  z  —  0  interface,  is 


,  _  B'(x,0)  Bi(x,0) 

•I n  —  ~  :  —  Jn  v 


■  l, 


Zs(*)  '  ^50 

Substituting  Equations  (2.20)  and  (A.l)  into  Equation  (A.  12)  yields 


(A. 12) 


A  = y-  r°°A  ‘,0rM  +  y-  r°°  \Am  +  simj  *i0rd0-  (-1.0) 

Z7T  J  oo  6 S ()  ^  -  OO  1 


Enforcing  the  fourth  condition  by  substituting  Equations  (A. 2)  and  (A. 13)  into 
( A.l  1 )  yields 


-  [ a'm 3)  -  «im |  +  t22-  |  Aim  -  elm | 

Kw>  +  fl»w)]-  ca.m) 

Equations  (A. 5)  -  (A. 9)  and  (A.l  1)  form  a  set  of  six  equations  that  can  be 
solved  for  the  six  unknown  functions,  '(/?),  .4,' ,(/?),  /?,\ (/i? ) ,  A^(^).  /?~(/?),  and 
Z?rj(/5).  Of  primary  interest  is  the  expression  for  the  electric  field  at  the  2  0 

interface 

Ei(*,0)  -  /  "J  I  .l!,(/3)  I  Dim  I  <M,T  <W-  (.1.15) 
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Noting  that  regions  k  =  0  and  k  —  3  are  the  same  media  and  solving  Equations 
(A. 5)  -  (A. 9)  and  (A. 14)  for  A\{f3)  and  Equation  (A. 15)  becomes 


i(M)  =  E”  F(0)j„  e^rf/3 

Z7T  J  —  OC 


(A. 16) 


where, 


F(0)  = 


/M7o  sinh(717’j)  +  /x07l  cosh^Tj) 


C>]70  b  JijjiJ fl07 1  sinh(7iTi)  +  [C#X(,7l  +  cosh(71T1) 


(A. 17) 


^  _  _j_  72  A270  cosh(7272)  A  /<-p72  sinh(72T2) 

■^50  2  [^270  sinh(72T2)  7  /io72  cosh(72T2) 


( A  A8) 


Now  consider  the  geometry  sketched  in  Figure  20(a).  The  geometry  is  a  sheet 
impedance  at  2  --  0  and  bordered  by  free  space  (/iy,ey)  for  z  >  0  and  material 
medium  (/ii,ej)  for  z  <  0.  This  geometry  can  be  considered  a  limiting  case  of 
the  two  layered  geometry  of  Figure  19  where  T\  — »  0  and  X2  — ♦  oo.  Since  region 
k  -  2  of  Figure  19  becomes  region  k  —  1  of  Figure  20(a)  then  the  k  —  2  subscripts 
will  become  k  —  1  subscripts.  Performing  these  limiting  operations  on  Equations 


(A. 17)  and  (A.18)  it  is  obtained  that 


FW  -  X71E 


1 _ l  )u  i  ;  l 

so 


(AA9) 


for  the  geometry  of  Figure  20(a).  If  the  sheet  impedance  is  entirely  in  free  space, 
as  in  Figure  20(b),  then  in  Equation  (A. 19),  — *  /ip  and  71  — >  79.  For  this  case, 


the  function  F{(3 )  is  given  by 


F(0) 


1 


'SO 


(A.20) 


This  ends  the  discussion  on  the  calculation  of  F(f3). 


61 


z 


(A6  0  ,  e  o) 


region:  k  =  0  f  Zs(x) 


region:  k  =  1 

(b) 


Figure  20:  (a)  Sketch  of  sheet  impedance  on  a  material  half  space,  (b)  Sketch  of 

sheet  impedance  in  free  space. 


APPENDIX  B 


THE  DETERMINATION  OF  FAR  ZONE  FIELDS  FOR  SEVERAL 
MULTILAYERED  SLAB  GEOMETRIES 


This  appendix  finds  the  far  zone  electric  field  of  current  expansion  function 
Jm  radiating  in  the  presence  of  the  multilayered  slab.  The  electric  field  expression 
of  Equation  (2.22)  for  region  k  —  0  is  evaluated  asymptotically  for  p  — >  o o  using 
the  method  of  stationary  phase.  The  result  is  given  for  the  three  geometries  illus¬ 
trated  in  Figures  19  and  20.  The  field  is  given  only  for  region  k  =  0  above  the 
multilayered  slab,  but  by  a  similar  approach,  the  field  below  the  multilayered  slab 
can  be  obtained. 


The  first  order  asymptotic  approximation  to  the  integral  [18] 

rb 


/(A)  =  f  f{t)e>X+Mdt  [BA) 

Ja 

for  A  — >  oo  and  t//(c)  =  0  with  a  <  c  <  b  and  c  not  too  close  to  a  and  b,  is 

(0.2) 


Employing  Equations  (2.22)  and  (A. 4),  the  field  for  region  k  —  0  above  the 
multilayered  slab  is 

£&(*,*)=  /h”  eith iff.  (»■:») 

J  —  oo 


To  get  this  expression  in  the  form  of  Equation  (B.l)  the  following  substitutions 
are  made: 


p  sin  <f) 


I2  ~  k?t  -  j  Jkl  f32 


yielding 


f  +  oc  ..  Jp  pros  <f>-  \  kl-(32  s\n4> 

eiM)=  /  A'!nw),-  I  v  I,/ 

J  —  oo 


Using  Equation  (B.2)  with 


<  -  f3 


m  -  w) 

— >  (3  cos  cj)  — 


ftsp  —  ~~  kf)  cos  0 


where  /3sp  is  the  stationary  phase  point,  it  is  obtained  that 


(HA) 


H,n(P,d>)  -  PmW 


(!<■  5) 


where 


Pm{4> )  =  ^(7n(-^0  cos  <f>)  J'l-Kky)  sin</>e;4  . 


(B.  6) 


Note  that  the  spectral  function  A^n(/3)  is  the  only  term  in  the  expression  for  E\ J, 
that  depends  on  the  multilayered  slab  geometry.  This  makes  it  straightforward 
to  find  the  far  zone  field  of  a  basis  function  radiating  in  the  presence  of  a  new 
multilayered  slab  geometry.  The  corresponding  to  the  new  geometry  must 

be  determined  and  then  inserted  into  Equation  (B.fi). 


fid 


The  function  A^n{(3)  can  be  found  by  enforcing  the  conditions  on  the  fields 


given  in  Section  2.4  and  following  an  approach  similar  to  that  outlined  in  Appendix 


A.  The  result  of  this  operation  is  that  for  the  two  layered  geometry  of  Figure  19 

Aim  = 

where 


1  / 

7* 

/X97ie7°Tl 

,CM170  +  J^P07t 

si«h(7iXi)  + 

Cp-071  +  ]^[Ml70 

cosh(7|Xi) 

c  =  -J-  +  -72 


A*270  cos  11(72X2)  +  /xo 72  sinh(72X2) 


(B.7) 

(B.  8) 


Zso  jut* 2  L/W270  sinh(72X2)  +  /z«72  cosh(72X2)  J  ' 

The  geometry  of  Figure  20(a)  can  be  considered  the  limiting  case  of  the  two 

layered  geometry  of  Figure  19  where  Xj  — >  0  and  X2  — >  00.  Since  region  k  =  2 

of  Figure  19  becomes  region  k  —  1  of  Figure  20(a)  then  the  k  —  2  subscripts  will 


become  k  =  1  subscripts.  Performing  these  limiting  operations  on  Equations  (B.7) 
and  (B.8)  it  is  obtained  that 

Aim  =  <B!» 

Zso  +  Wi 

for  the  geometry  of  Figure  20(a).  If  the  sheet  impedance  is  entirely  in  free  space, 


as  in  Figure  20(b),  then  in  Equation  (B.9),  /x  1  — *  /zo  and  71  — >  79.  For  tliis  case, 
the  function  4^(/3)  is  given  by 

Aim=  (B'I0) 

^50  ^  0 

To  obtain  a  closed  form  expression  for  the  far  zone  electric  field  above  the  mul¬ 
tilayered  slab,  the  function  x{(3)  is  evaluated  at  0  =  /3sp  =  —  kycoscf).  A^t(/9Sp) 
is  substituted  into  Equation  (B.6),  which  is  then  substituted  into  Equation  (B.5) 
yielding  a  closed  form  expression  for  Ern,  the  far  zone  field  of  current  expansion 


function  Jm.  Em  is  then  used  in  the  evaluation  of  the  voltage  vector  elements  of 


Equation  (2.36)  and  for  finding  the  total  far  zone  scattered  field.  This  ends  the 
discussion  on  evaluating  Em. 
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